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Abstract 

We reexamine the range of validity of finite-size scaling in the ip 4 
lattice model and the ip 4 field theory below four dimensions. We show 
that general renormalization-group arguments based on the renormal- 
izability of the ip 4 theory do not rule out the possibility of a violation 
of finite-size scaling due to a finite lattice constant and a finite cut- 
off. For a confined geometry of linear size L with periodic boundary 
conditions we analyze the approach towards bulk critical behavior as 
L — > oo at fixed £ for T > T c where £ is the bulk correlation length. We 
show that for this analysis ordinary renormalized perturbation theory 
is sufficient. On the basis of one- loop results and of exact results in 
the spherical limit we find that finite-size scaling is violated for both 
the ip 4 lattice model and the (p 4 field theory in the region L S> £. The 
non-scaling effects in the field theory and in the lattice model differ 
significantly from each other. 
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1 Introduction 



One of the fundamental achievements of the renormalization-group (RG) 
theory of critical phenomena is the elucidation and proof of universality and 
scaling near critical points [1-3]. These predictions have been shown to be 
asymptotically exact sufficiently close to criticality of infinitely large systems. 
For finite or partially finite systems, the field-theoretic version of RG the- 
ory has also provided an apparently exact prediction of universal finite-size 
scaling for systems with periodic boundary conditions || , in accord with phe- 
nomenological considerations || and with numerous analytical and numerical 
studies in statistical and elementary particle physics in the past decades [3, 
6-12]. Thus the validity of finite-size scaling appears to be well established. 

Consider, for example, the susceptibility x{t, L) of a ferromagnetic system for 
t = (T — T c )/T c > in a <i-dimensional finite geometry with a characteristic 
size L. For large L and small t the property of finite-size scaling means that 
X has the asymptotic form 

X (t,L)=x(t,oo)f{L/S) (1) 

where x(tj °°) = ^x^ -7 * s the bulk susceptibility and £ = £ot~ u is the bulk 
correlation length. For a given geometry and periodic boundary conditions 
the scaling function f(x) was found |4[] to be universal for d < 4 which implied 
that the relative deviation from bulk critical behavior 

s x(t,co) - X (t,L) = 

is universal as well in the entire range < L/£ < oo with g(oo) = 
1 — /(oo) = 0. This result, if correct, agrees with the scaling hypothesis p] 
which implies that the approach towards bulk critical behavior — > oo 
at fixed £ < oo above T c ) can be embodied in the scaling function g(L/^) 
|J. Universality in this context means that the shape of the scaling function 
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g(x) depends on the geometry and on the boundary conditions but does not 
depend on any nonuniversal parameter, in particular not on the lattice con- 
stant a of lattice models or on the cutoff A of field theories. As a consequence 
it was generally accepted that finite-size scaling functions such as g(x) can 
be calculated on the basis of field theories in the limit A — ► oo [3,4,7-50]. 

Brezin's RG analysis fij] started from the tp 4 lattice model with a finite lattice 
spacing a. The RG arguments, however, were presented within the renormal- 
ized theory after the limit a —>■ was taken. This limit is usually considered 
in studies of renormalized field theory of bulk systems for d < 4 [3, 13, 51] 
where cutoff and lattice effects are known to yield only subleading corrections 
to the leading critical temperature dependence. The asymptotic unimpor- 
tance of cutoff and lattice effects also for confined systems appeared to be a 
plausible assumption that was not questioned in Ref. Q] but is checked in 
the present paper and is found to be invalid. 



Very recently we have shown |52[ that this latter assumption is not generally 
justified in the 0(n) symmetric if 4 field theory of confined systems with pe- 
riodic boundary conditions. Specifically it was shown in the large-n limit for 
2 < d < 4 that a finite cutoff A implies a violation of finite-size scaling in the 
region L/£ ^> 1 above T c , with a non-exponential and nonuniversal approach 
A% oc (AL) -2 towards zero, even arbitrarily close to T c . This behavior was 
traced back to the (Vip) 2 term in the field-theoretic tp 4 Hamiltonian which 
approximates the more general interaction — <Pj) 2 of the (p 4 lattice model. 
In the latter model an exponential size dependence in the region L/£ ^> 1 was 
found |j2| , in accord with previous results for exactly solvable model systems 
[5,6,14, 53-57]. The possibility of a violation of finite-size scaling in the ip 4 
lattice model at finite lattice spacing, however, was not yet analyzed in our 
recent work f52 |. Thus the important question remained open whether the 



violation of finite-size scaling found in the (p 4 field theory [[5^] is an artifact 
of the field-theoretic continuum approximation or whether finite-size scaling 
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breaks down more generally for L/£ ^> 1 in confined lattice systems with a 
finite lattice constant. 

It is the purpose of the present paper to take up this problem for the O(n) 
symmetric <£> 4 lattice model at finite lattice constant a in the context of a 
detailed RG analysis, without taking the limit a —>■ 0. We assume renormal- 
izability in terms of bulk renormalizations and thus work for dimensionality 
d below the upper critical dimension which is 4 in our case. Thus this may 
become relevant to real three-dimensional systems. We shall show that the 
renormalizability of the ip A model in a confined geometry implies the asymp- 
totic (L — > oo,£ — > oo) validity of finite-size scaling for d < 4 at any fixed 
finite ratio L/£ < oo, in agreement with the proof of Brezin [f|, but does 
not rule out a violation of finite-size scaling in the limit L/£ — > oo at finite 
a/£ > 0. On the basis of one- loop results for general n and of exact results 
in the large-n limit we indeed find such a violation: instead of (2) the more 
general form 

A X = g(L/0 [1 + R(L/£, a/0] (3) 

must be considered where g(x) is universal but where the nonuniversal func- 
tion R contains a nontrivial dependence on the lattice constant a. Although 
R vanishes for fixed finite ratio L/£ < oo in the asymptotic region, 

lim R(L/£,d/0 = R{L/{,0) = 0, LH fixed, L/f < oo , 

(4) 

it exhibits a singular behavior in approaching the bulk limit L/£ — > oo at 
any fixed a/£ > 0, 

lirn^ R(x,d/^,) = oo , a/£ fixed, o/£ > . (5) 

This implies that for sufficiently large L/^ the leading size dependence 

A X ~ g{L/t)R{L/i,~a/t) , » 1 , a/^ > 0, (6) 

is nonuniversal and violates finite-size scaling. We emphasize that this vi- 
olation is not a subleading non-asymptotic property but occurs in leading 
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order at any finite £ < oo, even arbitrarily close to T c where "corrections 
to scaling" or corrections of bulk properties due to a finite lattice constant 
are completely negligible. (Here and in the following the symbol ~ means 
asymptotic behavior including the amplitude, i.e., A% ~ G(x) for x ^> 1 
means lim^oo Ax/G(x) = 1.) 

Our results imply that the property of finite-size scaling (for confined systems 
with periodic boundary conditions) which was previously believed to be exact 
for d < 4 in the asymptotic (L 3> a, £ ^> a) region of the L^ 1 — plane 
(Fig. 1) is not valid in a small but important part of this region (below the 
dashed line in Fig. 1). The nonuniversal function R is negligible at T = T c for 
sufficiently large L but it increases as L/£ increases above T c at fixed a/£ > 0. 
The approach to the bulk limit (arrow in Fig. 1) corresponds to a crossover 
from the scaling region to a non-scaling region where nonuniversal effects 
due to the finite lattice constant dominate the finite-size deviations from bulk 
behavior. The location of the (smooth) crossover region may be characterized 
by the line along which R(L/£,a/£) ~ 1, i.e., where the scaling and non- 
scaling contributions to Ax are equally large. This requirement defines the 
dashed line in Fig. 1. A similar line should exist below T c . Explicit results 
for R will be given in Sections 4 and 5 for the ip 4 lattice model for 2 < d < 4. 



Essential features of these results will remain valid even for d > 4 |58j such 
that the finite-size scaling form of Privman and Fisher |59|, Q will be violated 
for L ^> £. In particular, the lowest-mode approach of Brezin and Zinn- 
Justin |16[ and the phenomenological single-variable scaling form of Binder 



et al. |6T|, |62[ fail qualitatively for L ^> £ ^> a where these theories predict a 
universal power-law behavior Ax oc L~ d above four dimensions, rather than a 



non-universal exponential behavior [52, |58|| Ax oc e cL as derived in Sections 



4 and 5 of the present paper. Such striking structural differences [^j between 



the lowest-mode approximation and the effects of the higher modes cannot 



be regarded only as "corrections" [64 . 



4 



For comparison we also calculate the function Rfi e id(L/£,A£) for the field- 
theoretic (f A model at finite cutoff A (with periodic boundary conditions). 
For d < 4 we find a violation of finite-size scaling [^] due to a divergence 
of Rfi e id(L/£,A£) in the limit L/£ — > oo at fixed A£ < oo, analogous to (5) 
for the lattice model. The form of R field of the 99 4 field theory, however, 
differs significantly from that of R for the (p A lattice model. Even the sign 
of R f^id < is different from that of R > 0. Thus the y? 4 field theory based 
on the standard Landau- Ginzburg- Wilson continuum Hamiltonian does not 
predict the correct structure of the leading finite-size deviation from bulk 
critical behavior of lattice systems at any T > T c (and presumably at any 
T < T c ) for d < 4. We show that this statement remains valid also for 
d — 4 which may be relevant to elementary particle physics flO], [l^, [30], to 



disordered systems |65|] , and more generally to systems at their upper critical 
dimension . For d > 4 the failure of the continuum approximation is even 
more severe as it pertains to the entire L^ 1 — plane [58,63,66-68]. 

From a purely quantitative point of view, the non-scaling behavior of \ is 
a small effect that occurs predominantly in a region where the total finite- 
size contributions are exponentially small (for periodic boundary conditions). 
From a more fundamental point of view, however, the violation of finite-size 
scaling below four dimensions is a matter of principle, regardless how small 
this effect might be. In particular, our RG analysis for the simplest case 
of 'periodic boundary conditions raises considerable doubt about the validity 
of finite-size scaling in the more complicated cases of non-periodic boundary 
conditions where additional renormalizations and nonuniversal length scales 
come into play. They imply that general RG arguments are less compelling 
since additional assumptions would be needed. Furthermore, possible non- 
scaling effects for non-periodic boundary conditions may no longer be expo- 
nentially small. This may open up the prospect of resolving the longstanding 
and recent problems concerning the interpretation JTTJ] of experimental data 
for confined 4 He near the superfluid transition [69-80]. Work in this direction 
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is in progress. 



Our paper is organized as follows. In Section 2 we describe the renormaliza- 
tion scheme for the bulk (p 4 lattice model at finite lattice spacing. In Section 
3 we extend this scheme to the confined system and show that the RG argu- 
ments do not rule out a violation of finite-sized scaling in the limit L ^> £. 
In Section 4 we calculate A% for d < 4 in one-loop order for general n. The 
exact result for this quantity in the large-n limit is derived in Section 5. A 
summary and further discussion of our results is given in Section 6. 



2 Lattice Model: Bulk Properties at Finite 
Lattice Constant 

In this Section we introduce our notation and define the renormalization of 
bulk quantities of the if 4 lattice model at finite lattice spacing above T c . 
This will serve as the framework for the renormalization-group analysis of 
finite-size effects at finite lattice spacing in the subsequent Sections. A de- 
tailed formulation of the theory at finite lattice spacing is indispensable for 
clearly distinguishing lattice effects of the finite system from ordinary non- 



asymptotic Wegner |82fl corrections to scaling. 



We consider a ip 4 lattice Hamiltonian H for the variables ipi on the lattice 
points Xj of a simple-cubic lattice in a cube with volume V = L d and with 
periodic boundary conditions. (Generalizations to different geometries will 
be considered in the subsequent Sections.) We assume the statistical weight 
oc e~ H with 

H = ~ ad (e [tSV? + «o(tf ) 2 1 + E - (7) 

where a is the lattice constant. The variables (pi have n components (pi a with 
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a = 1, 2, • • • , n which vary in the range — oo < <p ia < oo. The couplings 
Jij are dimensionless quantities whereas the variables <fi have the dimension 
[a( 2_d )/ 2 ] and a has the dimension of a length. 



2.1 Unrenormalized Theory 

The renormalizat ion-group treatment of this model in the bulk limit V — > oo 
is well known which is usually carried out in the limit of zero lattice spacing 
a — > or in the limit of infinite cutoff A — ► oo in the continuum version (see 
(54) below) @, [U| [51]]. Here we shall formulate the renormalization of bulk 



quantities at finite a. Similar to the previous formulation of the bulk theory 
at fixed d < 4 |81] we shall express the bare theory in terms of the bulk 



correlation length £ before turning to the renormalized theory. 
We start from the bulk two-point vertex function |T3 , ^T|] 



T^(k,r ,u ,~a,d) = Xb(k)- 1 (8) 
where Xf>(k) is the bulk susceptibility at finite wavevector k above T c 

~2d 

Xb{k) = l/^ < > e****-**) . (9 ) 



It serves to define the bulk correlation length £ above T c according to 

e = ^(o^Mk)]- 1 | k ^ o • (io) 

We shall also consider the four-point vertex function of the bulk theory 



13, pB at vanishing external wavenumber. 



The parameter r in H is taken to be a linear function of the reduced tem- 
perature 

t = (T-T c )/T c , r = r 0c + a t , (11) 
7 



with ao > 0. The critical value r 0c of r is determined by Xfc(O) 1 = 0, i.e. 

T^(0,r 0c ,u ,h,d) = , (12) 
which provides an implicit definition of the function 

r 0c = r 0c (u ,a,d) (13) 

at finite lattice spacing a. (Note that ro c (wo, d, d) does not have an expansion 
in integer powers of u , unlike the vertex functions T^ N \r ,u ,d,d) |8l| .) 

Instead of r we shall substitute 

r = r - r Qc + r 0c (u , a, d) (14) 

into and consider (at k = 0) as well as the correlation length £ as 
functions of r — r 0c , Mo, d, d, 

r (A ° = rW(r -r 0c ,t 1 o,M), (15) 

f = £ ( r o - r 0c ,u ,a,d) . (16) 

Since £ is a monotonic function of r — r 0c , equation (16) can be inverted to 
define r — r 0c as a function of £, 

r -r 0c = h(£,u ,a,d) . (17) 

Finally we may substitute (17) into (15) which leads to bare vertex functions 
at finite a in terms of £, 

f^(e,«o,a,d) = rW(/ l (u,M),«o,M) ■ (18) 



These definitions are parallel to those at infinite cutoff in Ref. [EI]. (In 



particular, the functions have an expansion in integer powers of Uq.) 
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We illustrate these definitions by the one-loop results 



h(S,u Q ,a,d) = J„r l + 4(n + 2)!i J k (J k + J„r ) +O(u 2 



-1 



f( 2 )(k,e,u ,5,ci) = j k + j r 2 + o( M , 



r (4) (^,M ,a,rf) 
where 

with 



24m U - 4(n + 8)u / J k + J C 2 + O{u 2 ) 



k = [J(o) - j(k)j 



j(k) = x; j^-e-*-^-^-) 



(22) 



(23) 



We assume a finite-range pair interaction such that its Fourier transform 
has the small k behavior 



J k = J k 2 + 0{k\k] 



with a finite constant 



J = -(~a/L) d ^(Vo^Xi-Xj-) 2 >0 . 



(24) 



(25) 



The dependence of the quantities (19) - (21) on a comes from the integration 
limit of the bulk integral 

d k , _ _ N 

w (26 » 

K 

with \kj\ < 7r/a , j = 1, 2, • • • , d. Because of the super- renormalizability of 
the (f 4 theory || [13|, 5lJ the bare functions remain finite in the limit 
a — > at fixed £ and u for c? < 4. 
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2.2 Renormalization at Finite Lattice Constant 



As is well known, the perturbative results (19) - (21) of the bare theory do not 
provide a correct description in the critical region £ 3> a for d < 4. This prob- 
lem is circumvented by turning to the renormalized theory which provides 
a mapping from the critical to the non-critical region where perturbation 
theory is applicable [see (44) below]. We start from the bare iV-point vertex 
functions as functions of the correlation length £ where £ is considered 
to be a given quantity. The explicit determination of £ as a function of the 



reduced temperature t at finite a is a separate issue |8T| that is postponed 
to Appendix A. 

Since £ does not require a renormalization it suffices to introduce only two 
renormalization factors Z v and Z u to define multiplicatively renormalizable 
vertex functions fjj . We define the renormalized variable 

V? = Z; V Vi (27) 

and the renormalized coupling 

u = J 2 Z~ l u (28) 

with e = 4 — d. The reference length /i -1 is arbitrary. (It can be conveniently 
chosen as /i -1 = £ where £ is the amplitude of the asymptotic bulk corre- 
lation length as specified in (A. 17) of Appendix A.) The definitions (27) and 

~ (2) 

(28) lead to the following renormalized vertex function T R (at finite k) and 
f^P(at vanishing wavenumber) 

r%\k,£,u,ti,a,d) = Z^ f (2) (k, e, // Jl Z u u, a, d) , (29) 
fg ) (£,u,^,a,d) = Zlt^{^^JlZ u u,~a,d) . (30) 
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The Z factors Z v and Z u can be determined for d < 4 by standard renor- 
malization conditions at £ = 



dk 



d 

fg^k,// -1 ,^,//,^,^) 



k=0 



Jo 



r^(/x 1 ,u,n,a,d) = 24 j£ fi e u . 



(31) 



(32) 



The Z factors are finite for d < 4 if a > and remain finite for d < 4 if 
a — > at fixed u and /x. The following analysis is valid for d < 4 and d = 4 
since we keep the lattice spacing a finite. Substituting (20), (21), (24), (29) 
and (30) into (31) and (32) yields the Z factors in one-loop order 



Z v (u, fj,a,d) = 1 + 0{u 2 ) , 



(33) 



Z u (u,/id,d) = 1 + 4(n + 8) u I (an, d) + 0(u 2 ) (34) 

where 

I(an,d) = ^~ d J [n 2 + J k /J ]~ 2 • (35) 

k 

We derive renormalization-group equations (RGE) for f ^ by taking the 
derivative of (29) and (30) with respect to /i at fixed u , a and r — r 0c , i.e., 
at fixed £. This yields (at k = 0) 



N 



fidn + (3 u d u + — C 



with 



( v (u,fia,d) = (ndp InZ : ) , 



(36) 



(37) 



(3 u (u,fia,d) = (/i<9 M w) , 



(38) 
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where the index means differentiation at fixed parameters of the bare the- 
ory. The formal solution of (36) reads 

fSf } (£, u,im, a, d) = fW(£X^,M)exp^ J g(Oy (39) 

i 

where CyCO = (<p(u(£), ifia, d) and where u(i) is the solution of the flow 
equation 

= A,(«W,^,d) (40) 

with -u(l) = -u. In the present context the most convenient choice of the flow 
parameter £ is 

iiM = r 1 • (4i) 

We rewrite the renormalized vertex functions as 

f£° (e,«,^,a,d) = r 5 "/ w (/^^M) (42) 

where the amplitude functions are dimensionless and 

5jv = d-(d-2)N/2. (43) 

The renormalizability of the (p 4 lattice model for d < 4 guarantees that 
the limit a — > of at fixed and £ exists, i.e., that the function 
f^ N \^,u,0,d) is finite for finite /zf and u > for d < 4. From (39), (41) 
and (42) we obtain 

j\r r Apt 

fW(rt,u,vd,d) = /W(l,^),a/£,d)exp^ / ^(Oy (44) 

i 

which provides the mapping of the amplitude function f^ N \y,u,jj,a,d) from 
the critical region y ^> 1 (where perturbation theory breaks down) to the 
noncritical value y — 1 (where perturbation theory is applicable). The tem- 
perature dependence of the amplitude function f( N \l,u(l),a/£,d) in (44) is 
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affected by the finite lattice constant a not only through the explicit depen- 
dence on a/£ but also through the effective coupling u(£) that is determined 
by the a dependent RG flow equation (40). Furthermore the lattice constant 
enters the form of the temperature dependence of £(t) according to (A. 11) 
of Appendix A. 

2.3 Asymptotic Behavior 

Asymptotically (£ — > 0, £ — > oo) the effective coupling u(£) approaches the 
fixed point u* = u(0) as determined by 

= p u (u*,0,d) (45) 

which is independent of the lattice constant a and of the initial value u. For 
£ — > oo, equation (44) approaches the asymptotic form 

f^ N \^u,n~a,d) ~ AW /W(l,u*,0,d)(^) JV, ' /2 (46) 
with the critical exponent 

77 = -C„«0,d) (47) 
and the nonuniversal amplitude (which depends on -u and jj,a) 



A (N) = 




The amplitude function f( N \l, u, 0, d) is finite and nonsingular at u = u* > 
for d < 4 (compare the above statement after (43) regarding the renormal- 
izability of f^)- We see that the dependence on the lattice constant a 
has disappeared asymptotically (o/£ — > 0) in the amplitude function on the 
right-hand side of (46). Thus a enters the asymptotic bulk critical behavior 
of only via (beyond one-loop order) which is independent of £, and 
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via the amplitude £o of £ [see (A. 17)]. Therefore, taking the limit a — > 
in the renormalized quantity f ^ (as is usually done) is indeed justified in 
the asymptotic bulk theory since this limit does not change the asymptotic 
temperature dependence. For the asymptotic size dependence of the confined 
system, however, a corresponding property is not generally valid, as we shall 
see in the subsequent Sections. 

Application of the results to A^ = 2 yields, according to (8), the bare (phys- 
ical) bulk susceptibility at k = 

-l } df 

Xb = Z v (u,a/£ ,d) e[f (2) (l,u(i),~a/Z,d)]~ exp J (,(Oy ■ (49) 

e 

The asymptotic (£ — > 0, £ — > oo) critical behavior above T c is 

Xb = 4 = A+ t~~< . (50) 
The amplitudes depend on a according to 

4 = z v («,a/e ,rf)^ [A«/W(i,«*,o,d)]" 1 , (si) 

4 = ^4 , (52) 
where we have used \i = ^q 1 and the asymptotic form 

£ = 6 r 7 (53) 

with 7 = 1/(2 — 77). The dependence of £ 011 5 is given in (A. 17) of Appendix 
A. 
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2.4 Continuum Approximation 



For comparison we shall also consider the more standard version of the ip 4 
theory that is based on the continuum Hamiltonian 0, [U| |5l[ 



H 



d d x 



V 



\W + i(V^) 2 + Mo (^ 2 ) 2 



(54) 



for the n-component field y?(x). Here the (Vtp) 2 term approximates the 
(ipi — (fj) 2 term of the lattice Hamiltonian (7). The fluctuations of y(x) are 
confined to wavenumbers less than a finite cutoff A corresponding to %/a. 
The bulk susceptibility corresponding to (9) is now defined by 



Xfe(k) = lim f d d x < <p(x) <p(0) > 



-ikx 



(55) 



In the bulk limit, most expressions of the (p A lattice theory at finite a re- 
main applicable also to the ip 4 field theory at finite A after the replacements 
Jk — > k 2 , Jo — > 1 and a — > tt/A have been made. This implies that the 
asymptotic critical temperature dependence of the ip A field and lattice theory 
for bulk systems is identical (apart from different nonuniversal amplitudes). 
For confined systems, however, a corresponding statement regarding the size 
dependence is not generally valid as we shall see in the subsequent Sections. 



3 Renormalization Group and Finite-Size Scal- 
ing 

We are now in the position to discuss the size dependence of physical quan- 
tities within a renormalization-group treatment of the lattice model (7) with 
a finite volume V = L d . (For different geometries see below.) We focus our 
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analysis on the example of the susceptibility above T c 



X = 




<(Pi > 



(56) 



x(£,Uq,L, a,d) 



(57) 



Here we consider \ as a function of the bulk correlation length £ (rather than 
of ro) as explained in Section 2. 

The basic assumption in the following is that, for periodic boundary condi- 
tions, the ultraviolet divergences of x i n the limit a — > for the confined 
system are the same as those of the bulk susceptibility Xb- This plausible 
assumption is in accord with the argument of Brezin ||] regarding the de- 
composition of Fourier sums into bulk integrals (which carry the ultraviolet 
divergent part) and finite-size contributions (which are finite in the limit 
a — > 0). It is also in accord with one- loop results (Section 4) and with 
exact results in the large- n limit (Section 5). Although this assumption has 
far-reaching consequences regarding the validity of finite-size scaling for fixed 
finite L/£ we shall show that it does not rule out the possibility of a violation 
of finite-size scaling in the limit L>(. 

Our assumption implies that the renormalized susceptibility xr a t finite L 
and at finite a can be introduced for d < 4 as 



where Z^u, (j,a,d) and Z u (u, fj,a, d) are the bulk Z factors defined in Section 
2 and that xr remains finite in the limit a — > at fixed £,u,L and \x for 
d < 4. We note that there exists no justification a "priori, however, to actually 
perform this limit a — ► in the final results of the f A theory if they are to 
be compared with those of model systems on lattices with a finite lattice 



Xr(C,u,L,ii, a, d) 



Z v 1 x(£> // Jq Z ^ u, L, a, d) 



(58) 
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constant, e.g., with Monte Carlo data for Ising models. Therefore we shall 
keep a finite in the following analysis of the size dependence of xr an d of x- 

We derive a renormalization-group equation for xr by taking the derivative 
of (58) with respect to \i at fixed uq, a, L and r — r 0c , i.e., at fixed £. Since 
L is not renormalized M this yields 



Xr i€,u,L,ii,a,d) = (59) 



where p u (u,fj,a,d) and ( v (u, fAa,d) are defined by (37) and (38) of the bulk 
theory. The solution reads for d < 4 

j dl' 

XR{£,u,L,fi,a,d) = XR{£,u(£),L,£fi,a,d) exp J ^{l 1 ) — (60) 

e 

where i can be chosen arbitrarily in an exact theory. For the purpose of an 
application to perturbative results a natural choice is 

r 2 + l- 2 = fi 2 f . (6i) 

Although the RGE (59) has the same form as the bulk RGE (36) the simul- 
taneous appearance of the three lengths £, L and a in the arguments of xr m 
(59) and (60) complicates the situation and requires a careful consideration 
of different limiting cases. To answer the question about the possible rele- 
vance of the dependence on a we distinguish the following three cases (i) - 
(iii). 

(i) At T = T c or £ = oo and at finite L, we introduce a dimensionless function 
f x according to 

XR(oo,u,L,fi,d,d) = L 2 f x (oo, u, fxL, a/ L, d) . (62) 

In the large-L limit we obtain from (60) - (62) with i = n~ l L~ l 

-l ~ 



Xr{ 



oo, u, L, fi, a, d) ~ L 2 (fiL) 11 f x (oo,u* ,l,a/L,d) (63) 
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where r\ and are given in (47) and (48). The remaining dependence of 
/ x on a/L in (63) yields only a subleading correction to the leading power 
law ~ provided that the function / x (oo, u*, 1, 0, d) is finite. The latter 

property is valid for d < 4 (where u* > 0) provided that the f 4 theory is 
renormalizable at finite L and f x is non-singular at T c for finite fiL, i.e., 
provided that the limit (at fixed u and L) 

lim/ x (oo,ii, l,a/L,d) = f x (oo,u,l,0,d) (64) 

a— »0 

exists and is finite for u > at d < 4. 

(ii) For < £ < oo at finite L, i.e., at finite ratio < £/L < oo, we introduce 
the dimensionless amplitude function f x according to 

Xn(^u,L,fi,d,d) = £ 2 f x ((i£,u,fiL,a/L,d) . (65) 

In the asymptotic region £ ^> a, L ^> a corresponding to I < 1 we obtain 
from (60) 

«, L, /i, 5, d) ~ £ 2 f? IaW] ^ f x (£tf, u* , ifiL, a/L, d) . (66) 



Renormalizability of the ip A theory at finite L guarantees that the limit 

I™ a/L, d) = f x (f^,u,fiL,0,d) (67) 

a^U 

exists and is finite for finite arguments and d < 4. Therefore, taking the limit 
L — > oo in (66) at fixed finite ratio < L/£ < oo, i.e., at fixed values of 
and £//L, yields a finite amplitude function for d < 4 (w* > 0) 

f x m,u*,ef,L,0,d) = Y{L/i) (68) 

on the right-hand side of (66). Here we have used the fact that and ifiL 
depend only on L/£ according to (61). Thus the dependence on a/L in (66) 
represents only a subleading correction to the leading size-dependence (68) 
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provided that L/C is finite. This implies that both \ — Z v (u, /ia, d) xr an d 
Xr attain the finite-size scaling form 

XR (t, u, l, /i, a, d) ~ e^v-^i + e/L 2 r /2 [a®] ^ y{l/q m 

in the asymptotic region L>a and £ 3> a for any finite ratio < L/C < oo. 
This is in agreement with Brezin's conclusion Q who performed the limit 
a — > at the outset. 

(iii) There exist, however, significant paths in the asymptotic L -1 — £ -1 plane 
(Fig. 1) along which the ratio L/£ does not remain finite but diverges. These 
paths include the approach towards asymptotic bulk critical behavior at fixed 
a/£ > (arrow in Fig. 1). This case is not covered by the discussion of case 
(ii) above and was not considered in Brezin's analysis [|J. In this regime 
(L/£ S> 1, o/£ > 0) we make the choice 

t = ^r 1 (70) 

instead of (61). From (60), (66) and (70) we then obtain asymptotically 

Xr(£, u, L, fi, ~a, d) ~ e~>~ v [A (2) ] 1 f x (l, u*, L/C, a/L, d) . (71) 

In the bulk limit L — > oo at fixed £ < oo, equations (58) and (71) agree with 
(50) and (51) where 

/ x (l,«*,oo,0,d) = [/ (2) (l,M*,0,rf)] _1 = fZ . (72) 

Similar to the case (ii), the renormalizability of the (p 4 theory at finite L 
guarantees that in the limit a — > the function 

lim/ x (l,u*,L/£, a/L, d) = f x (l iU *,L/£,0,d) = f* b - h(L/C) (73) 

a— *u 

exists and is finite for finite arguments (with u* > for d < 4 dimensions). 
Thus, at first sight, the dependence of f x on a/L on the right-hand side 
of (71) appears to be a subleading correction that can be neglected in the 
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asymptotic region a/L <C 1, similar to the case (ii). A closer inspection 
shows, however, that this reasoning is not compelling in the present case (iii) 
where L/£ may become arbitrarily large at fixed finite £. 

On a formal level this is seen by rewriting the dimensionless function f x as 

/ x (l,«*,L/e,a/L,d) = F x (L/t,a/t) (74) 

where the dependence of F x on a appears in the form a/£ rather than a/L. 
Now there exists no argument why the dependence on a/( should be negli- 
gible in the limit L/£ 3> 1 at fixed a/£ > 0. More specifically, consider the 
decomposition 

F x (L/t,~a/0 = fb ~ fi(L/0 ~ / 2 (L/e,5/0 (75) 

where f b * and fi are independent of a as defined in (72) and (73). The last 
term f 2 contains the complete a dependence of f x and vanishes for a — > 0. 

Clearly only the term fi(L/£) in (75) is in agreement with finite-size scal- 
ing in contrast to the term / 2 (L/£, a/£). Thus the fundamental question 
arises whether the size dependence of f 2 at finite a/£ and for large L/£ is 
asymptotically negligible compared to that of fi. Only if the ratio 

RWtiK) = M f£)f < 76 > 

of equations (3)- (6) would vanish as L/^ — > oo at fixed < a/£ < oo there 
would be no violation of finite-size scaling. Renormalizability guarantees the 
existence of the function /i(L/£) but does not say anything about the mag- 
nitude and sign of the ratio i?(L/£,a/£) in the regime L/£ 3> 1. Although 
it is clear that the total finite-size deviation from asymptotic bulk critical 
behavior [compare (2)] 

Ax _ Xb-X = fb ~ f x {l,u*,L/£,a/L,d) (77) 
Xb fb 
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fi(L/Q 
fb 



1 + R(L/Z,a/Q 



(78) 



must approch zero as L/C, increases, renormalizability does not rule out the 
possibility that the relative contribution described by \R\ increases and be- 
comes large compared to 1 with increasing L at fixed £ < oo and a > 
in the region L >(. In fact it does not even rule out the possibility that 
R(L/£,a/£) diverges as L/£ — > oo at fixed < 5/£ < oo. If this is the case 
then /2 becomes dominant compared to f\ in (75) and finite-size scaling is 
violated in the lower part of the — plane close to the bulk limit (Fig. 

1)- 

We illustrate these considerations by a simple example: if f 2 (x,y) would 
be xy for small y and for general (arbitrarily large) x, then one should not 
dismiss f2{x,y) as a correction that is negligible (compared to fi(x)) for 
small y > since ft, can become large for x ^> 1/y. Actually we shall specify 
/ 2 (x, y) essentially as an exponential function of xy 2 [see equations (103) and 
(104)]. 

No general arguments but only explicit calculations can answer our question 
about the magnitude and sign of the nonuniversal quantity R. In Section 4 
we shall calculate R for general n in one-loop order. In Section 5 the exact 
form of R will be derived in the large-n limit. We shall indeed show that R 
diverges, 

lim R(x, a/£) = oo, a/£ fixed , a/£ > (79) 

at any finite a/£ for the (p A lattice model below four dimensions. Corre- 
sponding properties remain valid also in the cases d = 4 and d > 4 whose 



consequences will be studied elsewhere |SS 



The analysis of this Section can be repeated for the field-theoretic (p A model 
at finite cutoff A. The reasoning remains of course parallel to that given 
above and leads to the question about the magnitude and sign of the ratio 
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Rfieid{L/£,A£). We shall show in Sections 4 and 5 that Rfi e u differs fun- 
damentally from R and that in both cases finite-size scaling is violated for 

L/i » i. 

Finally we note that the analysis of this Section did not make explicit use 
of periodic boundary conditions except that the renormalizability in terms 
of bulk Z factors was assumed according to (58). Therefore our line of 
thoughts should remain applicable more generally to those cases where the 
bulk renormalizations suffice to renormalize the physical quantities of the 
confined system. Our analysis should also be extended to the important case 
where additional (surface) renormalizations come into play which, for con- 
fined systems, have been studied so far only in the continuum approximation 
(54) ( with surface terms) and only in the limit A — > oo. 



4 Perturbation Theory above T c 

As is well known a calculation of finite-size effects within the ip 4 theory 
including the size- dependence at T c requires a decomposition into modes 
where the lowest mode is separated and only the higher modes are treated 
perturbatively |I6|, [17|]. As noted recently |f)2| , however, this perturbative 
calculation does not correctly capture the exponential size dependence of the 
approach towards bulk critical behavior within the ip 4 lattice model. Covering 
the size dependence in the entire L^ 1 — plane, i.e., both at T = T c as well 
as for L/C, 3> 1 at fixed £ < oo, would require a non-perturbative treatment 



[ p^l within the mode expansion mentioned above. Such a treatment could 
be given on the basis of the order-parameter distribution function [ J7| that 
includes the higher modes in a non-perturbative way. 

Here we point out, however, that such a non-perturbative treatment can 
be avoided because the exponential size dependence of the ip 4 lattice model 
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above T c is not of a truly non-perturbative nature, unlike the effects due 
to the Goldstone modes below T c [f|7]] . We shall show that an ordinary 



renormalized perturbation approach is sufficient if one restricts oneself to the 
region L > £ above T c . This is just the region where the nonuniversal finite- 
size effects due to a finite lattice constant become significant. In this region 
a separation between the lowest mode and the higher modes is unnecessary 
since the lowest mode does not become dangerous in the bulk limit L — > oo at 
fixed £ < oo. It does not even become dangerous in the limit L —>■ oo, £ — > oo 
at fixed finite ratio L/£ > 0. Therefore we do not separate the lowest mode 
but instead shall present an ordinary perturbation approach, similar to bulk 
perturbation theory, where all modes are treated in the same way. Although 
this approach deteriorates with increasing £/L and breaks down in the region 
£ 3> L due to the dangerous lowest mode, it is well applicable for L > £. 



4.1 Lattice Model at Finite Lattice Constant 

We start out from the ip 4 lattice Hamiltonian (7) (for a d dimensional cube 
with V = L d and periodic boundary conditions) in the Fourier representation 



H 



V?kV?-k 



+ u L 3d (Wfu) (w'^-k-k'-k") > ( 80 ) 



kk'k" 



V?k = <>'> » ;k x r - ; (81) 

3 

where Jk is defined by (22) and (23) but now for the finite lattice. The 
summations in (80) run over discrete k vectors with components kj = 
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2nmj/L , rrij = 0, ± 1, ± 2, • • • , j = 1,2, ■■■,d in the range 
—7i /a < kj < 7i /a with a finite lattice spacing a. 

The standard one-loop expression for the (inverse) susceptibility (56) above 
T c reads for the finite system 

X' 1 = r + 4(n + 2) Mo L~ d E (r + J k ) -1 + O(ujj) (82) 

k 

and 

Xb 1 = r + A(n + 2)u J (r + J k ) _1 + O(ug) (83) 

k 

for the bulk system. Since J k is a periodic function of each component kj 
the one-loop sum in (82) satisfies the Poisson identity at finite a P, 83| 

L- d E (n> + ^k)- 1 = E / fa + e lknL (84) 

k n £ 

where k • n = J2j kj rij . The sum ^ n runs over integers rij, j = 
1,2, •■•,(! in the range — oo < rij < oo (whereas J2k an d Jk have finite 
cutoffs ±7r/S). The crucial quantity that contains all finite-size effects is the 
function 

D(r ,L,a) = L' d Efa + ik)" 1 - f (r + X)' 1 (85) 



E 

n^O k 



jo + 4 



-1 giknL 



(86) 



The sum X] m (85) includes the lowest-mode (k = 0) term L d r . It is only 

k 

this lowest-mode term that diverges for tq — > at finite L whereas the k ^ 
contributions remain finite in this limit. For finite r L 2 > 0, D(r Q , L, a) is of 
0{L 2 ~ d ) whereas for r L 2 3> 1 it is of 0(exp— L r^ 2 ) (see below). A simple 
rearrangement yields 



x(£,u ,L,a,d)~ 



J r 2 1 + 4(n 



2)u J 2 i 2 D 



+ 



0(t# 



(87) 
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with D = J D(J l; 2 ,L, a), i.e., 

D{t,L,a) = E(r 2 + 4/ Jo)- 1 I (r 2 + 4/ Jo)- 1 (88) 
k k 

where we have used Xb 1 = Jo + 0(uq) on the level of bare perturbation 
theory according to (8) and (20). 

The quantity D remains finite in the limit a — > for finite £ in arbitrary 
dimensions. This means that the ultraviolet divergence of x a ^ finite L in 
one-loop order is absorbed by the bulk correlation length £, in accord with 
the assumption in Section 3. For a — > and L/£ > 0, the function D can be 



represented as [16, 17 



D(£,L,0) = L 2 ~ d I{L 2 /e) (89) 

where 



d/2 



I(x) = (47T 2 )- 1 | tfe e""*/^ if(z) d - (~) 
o L 



(90) 



with 

oo 

K(z) = exp(-zm 2 ) . (91) 



m=— oo 



/(x) diverges for x — ► due to the lowest-mode term but is exponentially 
small for x ^> 1. For a > we decompose D as 

/J(£,L,a) = L 2 - d f D (L/£,~a/0 , (92) 

f D {L/£,a/0 = I{L 2 /i 2 ) + M(J/£ , 5/0 , (93) 

where M(L/£, a/£) contains the a dependence of D and vanishes for a — > 0. 
The explicit form of M(L/£, a/£) will be determined by equations (98), (101), 
(103) and (104) below. 

We note that D does not require a renormalization as it depends only on L, £ 
and a. Application of the RG analysis of the preceding Section to (87) yields 

X = Z v (u,na,d)f x (£jj^,u(£),£ijL,a/L,d)^ 2 exp / W)— (94) 
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where in one-loop order 



f x (tf,u,(iL,a/L,d) = J 1 {l-A(n + 2)u^r^/L) d - 2 f D (L/^a/O + 0(u 2 )} 

(95) 

For the application toL/£^>l at a/^ > we choose the flow parameter as 
£ = // £ . Then the finite-size deviation Ax from the bulk susceptibility 



Xb = Z^Ujfia, d)\ J 1 + O u(£) 



(96) 



becomes in one-loop order 

Xb — X „/„ , n\../«\/r-/r'i(i-2 r/r2 /t2\ 



Xb 
with 



A{n + 2)um/L) d - A I(L'/e) 1 + R{L/^a/i) 



+ O 



M(L/e,o/0 



(97) 
(98) 



These results are valid at finite a for both d < 4 and cZ = 4 dimensions and still 
contain all non-asymptotic contributions (Wegner corrections ||82|| ) within 
the ip 4 model. These contributions enter through u{£) as well as through 
the non-asymptotic form of £ as a function of t as determined by (A. 11) 
of Appendix A. Sufficiently close to T c such non-asymptotic contributions 
become negligible. By contrast, the nonuniversal term R(L/£,a/£) cannot 
be considered as a non- asymptotic contribution since it is nonnegligible at 
any T > T c for sufficiently large L as we shall see below. 

Neglecting the deviation of u(£) from the fixed point value u* we obtain for 
d < 4 

Xb - X 



Ax 



Xb 



g(L/0 1 + R(L/^d/0 



with the universal part 

g{L/Q = 4(n + 2) u* (£/ L) d ~ 2 I{L 2 / £ 2 ) + 0{u* 2 ) 



(99) 



(100) 



26 



For L/£ ^> 1 the function 7(L 2 /£ 2 ) becomes 

I(L 2 /e) = d^f-^iL/O^ 2 exp(-L/0 . 



(101) 



The nonuniversal part R(L/£,a/£) depends on the detailed form of J^. For 
simplicity we assume a simple-cubic lattice with nearest-neighbor coupling 
J, 



4J d 

= — ^2 i 1 - cos(akj)} 
a j=i 



(102) 



which implies Jo = 2J. A calculation parallel to that in Appendix A of 
Ref. |H yields for L > £ > a 



R(L/£,a/Z) = exp 



r(s/0 | 



with the function 



r(a/0 



24 



(103) 



(104) 



We see that, at any fixed a/£ > 0, i? diverges for L/£ — > oo. Thus the 
term R(L/£,a/£) is nonnegligible for sufficiently large L, even arbitrarily 
close to T c , unlike the Wegner corrections |82[ arising from the deviation 
of the effective coupling u(£) in (97) from its fixed point value u*. The 
resulting asymptotic size dependence of Ax for large L/£ and fixed a/£ > 
is determined by 



Xb 



g(L/Z)R(L/Z,a/Z) 



(105) 



4(n + 2) u* d (2tt 



,(l-d)/2 



(l-d)/2 



cxp 



1 

24' 



«/0 2 



L 
? 



exp(-L/0 

(106) 
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with u* > for d < 4. This behavior is nonuniversal and depends on the two 
ratios L/£ and a/^ in an essential way, rather than only onL/^. A dominant 
influence of o/£ exists in the non-scaling region where R > 1 corresponding 
to the region below the dashed line in the L~ l — £ _1 plane in Fig. 1. This 
line is determined by R = 1, i.e., 

a/L = [24U12]- 1 (d/^ 3 ■ (107) 

Non-negligible effects arising from R exist already above this line. 

The exponential part of (106) could be written in the form exp(—L/^ e ff) 
with £ e ff = £/[l — (a/£) 2 /24 + ...] which then would hide the violation of 
finite-size scaling [whereas the violation is quite explicit in the form of (103)- 
(106)]. But the true bulk correlation length £ (including all non-asymptotic 
bulk corrections) is already precisely defined by (10) and cannot be arbitrarily 
redefined here to comply with scaling. 

As a remarkable feature we note that, in one-loop order, the nonuniversal 
function R and the condition (107) are independent of the dimension d, of 
the number of components n and of the fixed point value u* of the four-point 
coupling. Therefore the same function R(L/£,d/£) causes a violation of the 
two- variable finite-size scaling form [59,60,63,66-68] above four dimensions as 
will be further discussed elsewhere J58| . 

These results can be easily generalized to a d dimensional system with a 
partially finite geometry that is confined in d dimensions (size L) and is 
infinite in d — d dimensions. This includes the cubic, film and cylindrical 
geometries as special cases d = d, d = 1 and d = d — 1, respectively. Instead 
of (82) we then have 

x -i = r + 4(n + 2)u L-*J2 [ ( r o + 4)" 1 + 0(v? Q ) (108) 
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where the <i-dimensional vector k = (q, p) has d components q = (qi, gj) 
and d — d components p = (pj +1 ), ■■■,Pd)- Equation (88) is modified accord- 
ingly. The integral representation of J(L 2 /£ 2 ), (90), is replaced by 



Kx) 



(4tt 



2\-l 



dz e - xz/ ^ 2 



K{z) 



7T 



i-d)/2 



(109) 



For L/£ ^> 1 this yields, instead of (101), 

/(L 2 /0 = d (27r)( 1 - d )/ 2 (L/0 (d ' 3)/2 exp(-L/£) . 



(110) 



The function R(L/£, a/£), however, remains unchanged, i.e., it is independent 
of the geometry (in one- loop order). Thus the leading finite-size deviation 
A% for large L/£ and fixed a/£ > for general d < d < 4 is given by 



X6 



4(n + 2)u*d (27r)( 1 - d >/ 2 (L/e) (1 " d)/2 



ex P ^(^/0 2 | exp(-L/0 • 

(Ill) 

This differs from (106) only by the replacement d — > d in the prefactor. 

We conclude that finite-size scaling is violated below four dimensions in the 
(p 4 lattice model with periodic boundary conditions above T c in one-loop order 
for general n in the region L/£ ^> 1 at any finite £ < oo even arbitrarily close 
to T c . Clearly higher-loop contributions cannot remedy this violation. 

Although one should trust the perturbative results of the (p A theory primarily 
for 2 < d < 4 one cannot exclude the possibility that an extrapolation to 
d = 2 yields sensible results, possibly for general n above T c and for n — 1 
below T c . This appears to be suggestive for the structure of our results (97), 
(106) and (111). This could be of particular relevance for the case n = 1 for 
which exact results of the two-dimensional Ising model in a confined geometry 
are available [5^, [84]. In all cases, contributions with an exponential size 
dependence were found above T c for L 3> £. It would be interesting to 



reanalyze these results p3, §4"| including all exponential and non-exponential 
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prefactors which so far have not been worked out explicitly (see, e.g., equation 
(6.8) of Ref. H) and to compare the structure of these results with that of 
our equations (106) and (111). Our results suggest that these prefactors may 
contain nonwaivers al contributions (such as our a dependent exponential 
factor) that cannot be neglected in the limit L/£ ^> 1 at fixed a/£ > 0. 



4.2 Field-Theoretic Model at Finite Cutoff 

The susceptibility above T c of the field-theoretic model (54) is defined as 

X = J d d x < p(xMO) > . (112) 

v 

The one-loop expression of the susceptibility for a cubic geometry, V = L d , 
reads [compare (87) and (88)] 



X(£, u , L, A, d)- 1 = r 2 [l + 4(n + 2)u £ 2 J D /ieM + O(u 2 )} (113) 
where now 

Dfieid(£, L, A) = L- d J2 (r 2 + k 2 )- 1 - J (r 2 + k 2 )- 1 . (114) 



k 



k 



Here the range of k is limited by the cutoff A according to \kj\ < A for the 
bulk integral / and —A < kj < A for the sum J2k- F° r A — > oo the function 

k 

D field becomes identical with D for a — » 0. Thus we decompose 



n — r 2 ~ d 

u field — L 



I(L 2 /e) + M field (AL, AO (115) 



with Mf ie id(oo, oo) = where the function I(L /£ ) is the same as for the 
lattice model for a — > 0, (90) and (91), but the cutoff dependent part M/j e zd 
of D fi e id differs fundamentally from the a dependent part M of D, as we 
shall see. 
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Application of the RG analysis of Section 3 to the field-theoretic model leads 
to 



Xb-X 
Xb 

with 



A(n + 2)u(£)(aL) d - 2 I(L 2 /e)[l + R field (L/^ AL)\ + 0[u(t)' 

(116) 



R field {L/^AL) 



M field (AL, A£) 



(117) 



/(Lve) 

Equation (116) is valid for d < 4. For AL ^> 1 and A( > 1 we have found 



M /ieW (AL,A£) 
where 

a (d) 



d a (d)(AL) d - 4 + O (AL) d -\exp(-A- 2 £ 



2t-2\ 



— I dxx e x 



i 

— / dy exp(-y 2 x) 

-1 



d-1 



(118) 



119) 



Together with 7(L 2 /£ 2 ), (101), this yields the large LA behavior of Rfi e i d at 
fixed A£ > 1 



R field {L/i,AL) 



i2nY d - 1 ^ 2 a (d)(L/O {3 ' d)/2 (AL) d - i e L/ «. (120) 



We see that, at fixed A£ ^> 1 , Rfi e i d diverges exponentially towards — oo 
for L/£ —>■ oo and Ax has the asymptotic size dependence in this limit (for 
cubic geometry and d < 4) 



Xb-X 
Xb 



g(L/Z)R fiM {L/Z,AL) 



-4(n + 2)u*d a (d)(A£) d - 2 (AL) 



;i21) 



(122) 



For the non-cubic geometries defined in Subsection 4.1 the corresponding 
results are obtained by substituting 7(L 2 /£ 2 ) in the form of (110) instead of 
(101) and by the replacement d — > d in the prefactor of (122). Rfi e u remains 
unchanged. 
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Like the result (106) for the lattice model, the behavior (122) is nonuniversal 
and violates finite-size scaling, as pointed out recently f52f . The non-scaling 
effect becomes significant in the region below the dashed line in Fig. 1 which, 



for the field-theoretic model, is determined by \Rfi e u\ = 1 |52[] , 

(2nf- d ^(L/0 {d - 3)/2 exp(-L/0 = a (d)(AL) d ~\ (123) 
for both cubic and non-cubic geometries. 

The structure of (122) differs fundamentally from that of (106) for the lattice 
model in three respects: (i) the size-dependence of (122) is non-exponential, 
(ii) the dependence on A is non-exponential, and (iii) the sign of (122) is 
negative. Since for ( > L we must have Xb — X > this sign implies the 
existence of a crossing point of the bulk curve for Xb and the curve for x a ^ 
some T > T c , in contrast to (106) where x does not cross the bulk curve for 
T >T C . We conclude that the y? 4 field theory does not correctly describe the 
leading finite-size deviations from the bulk critical behavior of lattice systems 
below four dimensions. This is true also at d > 4; for d = 4 this follows from 
(97), (103) and (116), (120). The case d > 4 will be discussed elsewhere 



In Section 3 we have employed renormalization conditions in order to define 
the renormalized theory at finite a and finite A for d < 4. For a calculation 
of the universal part g(L/£) of the finite-size effect, however, it is possible to 
employ a more convenient RG approach using dimensional regularization and 
minimal subtraction at fixed 2 < d < 4 ||8X|, |85[] , as has been done in the finite- 
size calculations of Refs. [39-41,46,47,50]. In this case the fixed point value 
u* in (106) and (122) is replaced by A^u^in where now u* min can be taken 
from the accurate Borel-resummed results of the minimally renormalized bulk 
theory [El], FJJ ^7|. For d = 3 the corresponding values are A% 1 = An and 
u* min = 0.0404, 0.0362, 0.0327 for n = 1,2,3, respectively |87|. In the large 
71,-limit (at fixed un) the exact fixed point value is u* min n = (4 — d)/4 for 
2 < d < 4. 
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5 Exact Results in the Large-n Limit 



In the following we perform the analysis of finite-size effects in the ip 4 model 
above T c in the large-n limit at finite lattice constant and finite cutoff for 
d < 4 without using the renormalization group. The exact results in this 
limit will confirm the perturbative RG results of the preceding Section, in 
particular the existence of a non-scaling region for both the field-theoretic 
and the lattice <£> 4 model in the range L ^> £. The existence of this region 



was overlooked in our recent work [56 



5.1 Lattice Model at Finite Lattice Constant for n 



oo 



We start from the exact result for the susceptibility x/ n = X P er component 
of the ip A lattice model above T c in the limit n —>■ oo at fixed uon in a cubic 



geometry as determined by the implicit equation [66 



X 



r + Au nL d ^(J k + x)~ l . 

k 



For d > 2 this can be rewritten as 
XT 1 = r -r 0c + 4M nD(x _1 ,^,a) - ^UqTix' 1 

where L,d) is defined by (85) and 

r 0c = -Au n / J k 2 



X 



-1\ 



-1 



(124) 



(125) 



(126) 



The bulk susceptibility Xb an d the bulk correlation length £ are determined 
by 



X b 



r - r 0c - Au nx b 1 / Jk (X + Xb 



(127) 
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£ 2 = JoXb 



(128) 



Although RG arguments will not be needed in the following we note that the 
ultraviolet (a — > 0) behavior of x > (125), is the same as that of Xb because 
the function D has no ultraviolet (a — > 0) divergence. This supports the 
assumption made in Sect. 3. In the following we keep a finite. 

For L ^> a and £ a, i.e., for small x~ l a 2 at finite a, the bulk integral in 
(125) yields for 2 < d < 4 

/ [4 (4 + r 1 )]" 1 = ^ 4, X e/2 e- 1 {l + O [(r'-aY 2 ]} (129) 



with e = 4 — d and 



A d = 2 2 - d 7r- (i/2 (rf-2)- 1 r(3-rf/2) 



For u nJ d ^ 2 A d e 1 x e ^ 2 > 1 this leads to 



X = Xb 



eAf i^b{j]l 2 tl\Wa) 



2/(2-<2) 



(130) 



(131) 



where the function D is defined by (88). For L>(we may replace Jo^X 1 ^ 2 
by Jl^xJ 2 — £ in IX Using the decomposition (93) we arrive at 



Ax 



Xb-X 
Xb 



= g(L/t) l + R(L/Z,a/Z) 



where 



g(L/0 = 2eA d \d-2)-\i/L) d ^I{L 2 /e) 



with I(L 2 /^) given by (101), and 



R(L/S,a/£) = exp 



- 1 



(132) 



(133) 



(134) 



with T(a/£) given by (104). The structure of g agrees with that of g, (100). 
(We note that the factor eA^ 1 in (133) can be interpreted as Au*n where u* 
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is the fixed point value in the large-n limit, see the last paragraph of Section 
4.) The function R turns out to be same as R, (103) and (104), which was 
to be expected because R is independent of n (in one- loop order). 

For non-cubic geometries defined in Section 4.1 the result is only modified 
by the replacement d —>■ d in the prefactor of the expression for 7(L 2 /£ 2 ) as 
given by (110), thus g reads explicitly 

g{L/i) = 2 d tt 1/2 [r(e/2)]- 1 (2e/L) (d - 1)/2 e- L/? . (135) 

Equations (132)- (135) prove that finite-size scaling is violated for 2 < d < 4 
in the large-n limit for L/£ ^> at fixed a/£ > above T c . This exact result 
supports the correctness of our conclusions in the preceding Section based 
on one- loop results. We note that the results (133)- (135) have finite limits 
for d — > 2 at fixed £. 



5.2 Field-Theoretic Model at Finite Cutoff for n 



oo 



The analysis of the susceptibility x = x/n in t ne field-theoretic <£> 4 model at 
finite A for n —>■ oo at fixed u^n is parallel to that in Section 5.1. Equations 
(124) - (130) remain valid after the replacements — > k 2 and J ~~ * 1- 
Instead of (131) we now obtain for L > A -1 and £ » A -1 (for 2 < d < 4) 

2/(2-d) 



X 



Xb 



1 + eA^e- 2 D field (^L,A) 



where is defined by (114) and (115). For L > ^ we arrive at 



Ax 



Xb -x 

Xb 



g(L/0 l+R field {L/£,AL) 



(136) 



(137) 



where g is identical with (133) or (135). Furthermore Rfieid is identical with 
R field, (120), as expected because Rfieid is independent of n (in one- loop 
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order). These exact results confirm the conclusions drawn in Section 4 about 
the violation of finite-size scaling in the ip A theory and about the failure of 
the continuum approximation for finite lattice systems in the region L ^> £ 
above T c . 



5.3 Spherical Model 

Since the ip A lattice model in the large-n limit and the spherical model []88|, 
are expected to yield asymptotically identical results (for the case of periodic 
boundary conditions) we compare our results with those by Barber and Fisher 
[54j1 and by Singh and Pathria |56| . Here we only study the approach of the 



susceptibility towards the bulk critical behavior above T c . 

Barber and Fisher (BF) considered the spherical model with a film geometry 
which should correspond to our result (132), (134) and (135) with d — 1. The 
result of BF for the dimensionless susceptibility xbf reads for fixed T > T c 
and L/a — > oo 

Xbf ~ Xbf + J^Bl{T){L/~a)^/ 2 eM-^d{T)L/d] (138) 

where Xbf ^ s the dimensionless bulk susceptibility and Jo = 2J. (Here we 
have corrected a misprint in (8.11) of Ref.[[54[] by replacing the exponent 
(3 - d)/A by (1 - d)/2, compare (8.9) and (8.10) of Ref.g].) Both functions 
B®(T) and r^(T) are expressed in terms of the dimensionless function $o(^) 
which is given by the dimensionless inverse bulk susceptibility according to 
&o(T) = Jq {Xbf)~ 1 = (<V£) 2 - Calculating the small-$o behavior of the 
derivative of the generalized Watson function H^($ ) as 

W' d {%) ~ -2- d n- d/2 r(e/2) $ e/2 (139) 

with e = 4 — d we have found 

Bf\T) ~ -2( d+1 )/ 2 7r 1 / 2 [r(e/2)]- 1 ^ d+3 ^ (140) 
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for T — > T c . This leads to 
Xbf ~ Xbf 

X B F 



2 vr 1 / 2 [r(e/2)]- 1 (2^/L)( d - 1 )/ 2 e- r ^ T ) i / s . (141) 



The non-exponential part agrees with that of our g(L/Q, (135), for d = 1. 
If we expand the exponent T^(T) = 2 arcsinh($,!/ 2 /2) to third order in ^\j 2 
and express it in terms of a/£ we obtain 



exp 



T d {T)L/~a 



exp 



(2L/a) arcsinh ( ^/C 



(142) 



exp 



L/C 



24 



(a/0 2 L/e + O (a/0 4 V£ 



(143) 



which also agrees with the exponential part of our solution (132), (134), (135) 
of the lattice model. The importance of the positive second term oc (a/£) 2 
in (143) for the leading finite-size deviations from bulk critical behavior (for 
L/C — ► oo at fixed a/£ > 0) was overlooked by BF [|54| (who considered the 
exponent in (138) only in the limit T — > T c at fixed Ljd < oo rather than 
L/a — > oo at fixed T > T c , see also equations (5.75)-(5.77) of Ref.|^]). 



The solution of Singh and Pathria (SP) for A% (see equation (19) of Ref. 
[|5f|) in the large-L limit at fixed T > T c agrees with the universal part g of 
our result, (135). Our non-universal contribution R, (134), is not contained 
in the solution of SP. 



6 Summary and Conclusions 



In the following we summarize and further comment on the results of this 
paper as follows. 
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We have studied the consequences of renormalizability (in terms of bulk 
renormalizations) of the ip 4 theory in a confined geometry with periodic 
boundary conditions below four dimensions. We have found that the conse- 
quences for the confined system are in contrast to those for the bulk system. 
While for the bulk system renormalizability implies that the leading critical 
temperature dependence is not significantly affected by a finite lattice con- 
stant a or by a finite cutoff A this is not generally the case for the leading size 
dependence of the confined system. Lattice and cutoff effects are asymptot- 
ically negligible as criticality is approached at fixed finite ratio L/£ but not 
in the limit L/£ 3> 1 above T c . In the latter case the leading finite-size effect 
on the susceptibility \ turns out to be nonuniversal, i.e., to depend explicitly 
on a and A [see equations (106) and (122)], and to violate finite-size scaling 
in the region L ^> £ of the L^ 1 — plane (Fig. 1). 

Although lattice and continuum models yield the same asymptotic critical 
behavior of bulk systems we have shown that this is not the case for confined 
systems. While lattice systems (with periodic boundary conditions) have 
an exponential size dependence of A% above T c for large a power-law 
behavior Ax oc (AL)~ 2 is obtained from the ip 4 field theory ||52|| . 

It is expected that part of our conclusions apply also to T < T c for the 
case n = 1. For n > 2 (and for periodic boundary conditions) the power- 
law behavior of the Goldstone modes is expected to govern the finite-size 
deviations from bulk critical behavior such that nonuniversal exponential 
terms become subleading . 

Our results for the leading finite-size deviations from bulk critical behavior 
have been derived for 2 < d < 4 dimensions but our conclusions may be 
applicable to d = 2 dimensions. A detailed reexamination of the exponential 
terms in the existing exact results for the two-dimensional Ising model |53|, p4| 
would be interesting. 
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Part of our results remain valid also at and above four dimensions. One of 
the consequences is that the universal two-variable finite-size scaling form 
for the y? 4 lattice model for d > 4 [59,60,63,66-68] is violated for L/£ ^> 1 
above T c for general n and below T c for n — 1. A further consequence is that 
the predictions of the lowest-mode approach [16) and of the phenomenological 



scaling theory implying a single- variable scaling form |61, fail qualitatively 
in the region £/L ^> 1 where these theories predict a universal power-law 
behavior Ax oc L~ d instead of a non-universal exponential size dependence 
Ax oc e~ cL . The latter can easily be incorporated in our theory above four 
dimensions |63| by extending our present perturbation approach of Section 4 



to d > 4 58 



From a quantitative point of view, our prediction of a violation of finite-size 
scaling is difficult to be tested by means of Monte Carlo simulations (e.g., 
for Ising models) because the non-scaling effect on x occurs predominantly 
in the region L>( where the total finite-size effects on x are exponentially 
small (for periodic boundary conditions). 

Our general arguments regarding the consequences of renormalizability (as 
presented in Section 3) are presumably not restricted to periodic boundary 
conditions but may be generalized to non-periodic boundary conditions. We 
consider this to be potentially important for applications to real systems 
where finite-size deviations from bulk critical behavior are not exponentially 
small. We do no longer see a stringent reason to believe that renormalizabil- 
ity implies the validity of finite-size scaling in the more complicated cases 
of non-periodic boundary conditions where additional renormalizations and 
nonuniversal length scales come into play. In particular for the important 
case of confined 4 He near the superfluid transition where the entire region 
L > £ and L < £ is perfectly well accessible to high-resolution experiments 
[69,70,72-77,79,80,90] we cannot exclude the existence of nonuniversal non- 
scaling effects in the (f 4 (lattice and field) theory with Dirichlet boundary 
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conditions. This could eventually lead to a natural explanation of longstand- 
ing and recent discrepancies between experimental data [69-80] and theoret- 
ical predictions [5,11,24,29,31,32,36,37,42,49] that were based on (seemingly 
plausible) assumptions which imply the validity of finite-size scaling. Also for 
the exploration of finite-size effects on transport properties in 4 He on earth 



[72,74,79,90] and under microgravity conditions [91] as well as on thermody- 



namic properties near ordinary critical points under microgravity conditions 



92], detailed knowledge on the effect of a finite atomic distance may turn 



out to be important. 
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Appendix A: Bulk Correlation Length at Fi- 
nite Lattice Constant 



In this Appendix we determine the bulk correlation length £ above T c as 
a function of t at finite lattice spacing for d < 4. The derivation is par- 



allel to that at infinite cutoff in Ref. We introduce the renormalized 

temperature variable 

r = Z~ x (r - r 0c ) = Z~ l a n t 



(A.l) 



where Z r is identical with the Z factor Z^i 



Z r that is needed to renormalize 
the bare vertex function n 1 - 2 ) (£, u , a, d) ||, [L3|, |5l|] for d < 4 

f £' 2) (£, it, /i, a, d) = Z r Z v (£, // J 2 Z u it, a, d) . (A.2) 



The one-loop expression for f ( 1>2 ) is 

f^ 2 ) (e,«o,a,d) = 1 - 4(n + 2)«o/ + 



+ 0( M 2 ). (A.3) 



The Z factor Z r (u,/j,a,d) for d < 4 is determined by the renormalization 
condition at £ = /i -1 



r 



(1.2) 
R 



(/i 1 ,u,fx,a,d) 



1 



which yields in one-loop order 

Z r (u,/ia,d) = 1 + 4(n + 2) ul(fid, d) + 0(u 2 ) 



(A.4) 



(A.5) 



where I(fia,d) is given in (35). Using (A.l) we rewrite the right-hand side 
of (17) in terms of renormalized quantities as 



r - r 0c 



h(£, jJ, e J Q Z u u, a, d) 



(A.6) 



Z r (u, fxa, d) /i 2 Q(/i£, u, fid, d) 



(A.7) 
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with the dimensionless amplitude function Q. Taking the derivative at fixed 
u , a and r — r 0c (i.e. fixed £) yields the RGE 

+ AA + (2-Cr)] Q{rf,u,iM,d) = , (A.8) 



( r (u,/ia,d) = (fi df, In Z r 1 ) . (A. 9) 

The formal solution is 

Q(j4,u,fi,a,d) = Q(l,«(*),a/£,d)exp y (C - 0(0] -y (A.10) 

1 

with £ = (/^O -1 where (* = ( r (u*,0, d) and £r(0 = Cr{u(i),£^d, d). 
Equations (A.l), (A. 7) and (A. 10) can be summarized as 

T dt> 

r = at = C 2 Q(l,u(e),~a/Z,d) exp J ( r (£') — (A.ll) 

i 

with 

a = Z r (u, fid, d) _1 ao > (A. 12) 

where the bare parameter ao is defined in (11). Equation (A.ll) determines 
t > as a function of £ for d < 4 at finite a. Inverting (A.ll) yields 
£(t) including non-asymptotic (Wegner |82]| ) corrections. The asymptotic 
(£ — > oo) form of the correlation length follows from (A.ll) as 

£ = tot-" (A. 13) 

with the critical exponent 

v = (2-C)- 1 (A.14) 

and the amplitude 



to = fi a 



Q*exp j [Q -Cr(t')) (A. 15) 
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where 



Q* = Q(l,u*,0,d) . (A.16) 

After the choice \i = the correlation- length amplitude £ is determined 
implicitly in terms of the bare parameter ao, the lattice spacing a and the 
renormalized coupling u by 

°r df 
e = Z r (u,~a/to,d)a 1 Q*exp J [£ - ( r (£')) — . (A.17) 

1 

We note that the functions h(£, Uo,a, d) and Q(/j^,u, fj,a,d) do not have an 
expansion in integer powers of u and u, respectively, beyond one-loop order 



The one-loop expression of Q(n£, u, /xa, d) and of Q* can be derived 
from (19), (34), (A. 5) and (A. 7). An integral representation of Q in terms of 
expandable functions can be derived at finite a along the lines of Section 4.1 
of Ref. 
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Figure Caption 



Fig. 1. Asymptotic L^ 1 — plane in units of the lattice constant a above 
T c (schematic plot) for the <£> 4 lattice model below four dimensions where L 
is the system size and £ is the bulk correlation length. Nonuniversal lattice 
effects become nonnegligible in the non-scaling region below the dashed line 
where finite-size scaling is violated. This crossover line is determined by 
equation (107). Well above this line the dependence on the lattice spacing a 
is negligible in equation (3). The arrow indicates an approach towards bulk 
critical behavior at constant < t <C 1 through the non-scaling region where 
equation (106) is valid. A corresponding non-scaling region exists also for the 
field-theoretic model at finite cutoff where the crossover line is determined 
by equation (123) and the non-scaling effect is described by equation (122), 
compare Fig.l of Ref. [52]. A corresponding crossover line should be added 
to Fig.l of Ref. [63]. 
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